What is claimed is: 



1. A method of model reduction and system identification of a dynamic system 
with multiple inputs, comprising: 

generating a plurality of statistically independent random numbers for use as 
input signals; and 

performing a singular-value-decomposition directly on a system response of 
the dynamic system due to a simultaneous excitation of the plurality of input 
signals. 

2. The method of Claim 1, further comprising sampling individual pulse 
responses for a first time step and a second time step. 

3. The method of Claim 1, further comprising sampling a system response y" for 
n = 0, 1, 2, . . . M. 

4. The method of Claim 1, further comprising defining Hankel-like matrices H c0 
and H c i as follows: 
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Hco = bio yf Q --• y^Tl 

= Clx 1 xr ... x jW - 3 i (25) 

Hci - Hi yfj - y^-H 

= CAbt 1 x 2 ... x* w - ! | (26) 

Hco ■ USV T 

= U^S^VS (27) 

5. The method of Claim 4, further comprising obtaining system matrices (A, B, 
C, D) by a least square approximation as follows: 

D = Y 0 (28) 

C a U.KSjf (29) 

B - S^U^Y 1 (30) 

A S^^USHdV^S^' 2 (31) 

6. The method of Claim 4, wherein (M- 1) and No 

7. The method of Claim 4, wherein a total number of input samples is equal to 
M+ 1 +2xN h 

8. The method of Claim 1, further comprising defining augmented Hjoi and H c n 
matrices as follows: 
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where; 
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Afx 1 xr ... it* 1 - 1 ] 



JC2 



(32) 



(33) 



10 



9. The method of Claim 8, wherein a total number of input samples is equal to 
M+l+K+(2+K)xNi. 

10. The method of Claim 1, wherein at least some of the input signals are filtered 
through a low-pass filter. 

11. The method of Claim 1, wherein the plurality of input signals includes 
applying multiple step inputs in a sequential manner. 

12. The method of Claim 1, wherein the plurality of input signals includes 

applying multiple pulse inputs in a sequential manner. 
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13. The method of Claim 1, further comprising performing a second order 
reduction based on the Frequency-Domain Karhunen-Loeve (FDKL) method 
to the SCI/ERA ROM using the plurality of input signals. 

14. The method of Claim 13, further comprising premultiplying the SCI/ERA 

T 

ROM matrices by * to yield a new reduced-order model as follows: 

p® 4L - AipF+BamP (S3) 

= CxpP + Du" 054) 

where 

A a s # r A* $5) 

B a e * T B (56) 

Ci = C* (ST) 

15. A method of model reduction and system identification of a dynamic system 
with multiple inputs, comprising: 

generating a plurality of statistically independent random numbers for use as 
input signals; and 

performing a singular-value-decomposition directly on a system response of 
the dynamic system due to a simultaneous excitation of the plurality of input 
signals; 

sampling individual pulse responses for a first time step and a second time 
step; 

defining H c0 and H c i matrices as follows: 
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SVDofHuQ yields 



Hco = bio Yea ■■■ y'aTh 

= C|k 1 ai 3 ... x^U (25) 

Hci = bit yfa — yci~ L !i 

= CAf* 1 x 3 ... x*- 1 ! (26) 

Hco = USV r 

= U*Sjf 2$ 2 V£ (27) 



; and 

obtaining system matrices (A, B, C, D) by a least square approximation as 
follows: 

D = Y° (28) 

C a (29) 

B - E^-U^Y 1 (30) 

A ^ S^^USH.xY^S^ 2 (31) 

16. The method of Claim 15, wherein at least some of the input signals are filtered 
through a low-pass filter. 

17. The method of Claim 15, wherein the plurality of input signals includes 
applying multiple step inputs in a sequential manner. 



18. The method of Claim 15, wherein the plurality of input signals includes 
15 applying multiple pulse inputs in a sequential manner. 
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19. 



The method of Claim 15, further comprising performing a second order 
reduction using the Frequency-Domain Karhunen-Loeve (FDKL) method to 
the SCI/ERA ROM using the plurality of input signals. 



20. A method of simulating a fluid flow, comprising: 

generating a plurality of statistically independent random numbers for use as 
input signals; and 

performing a singular-value-decomposition directly on a fluid response due to 
a simultaneous excitation of the plurality of input signals. 

21. The method of Claim 20, further comprising sampling individual pulse 
responses for first and second time steps. 

22. The method of Claim 20, further comprising defining H c0 and H c i matrices as 
follows: 



Hco s h r lo Yen -■- Jcd" 1 ! 



(25) 



(26) 



SftDafH^ yields 



USV T 




(27) 
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23. 



The method of Claim 22, further obtaining fluid system matrices (A, B, C, D) 
approximately as follows: 



(28) 
(29) 
(30) 
(3X) 



24. The method of Claim 22, further comprising denning augmented H c0 i and 
H c ii matrices as follows: 
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where 



JVj 



(32) 



(33) 



j-t 
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25. The method of Claim 20, wherein at least some of the input signals are at least 
one of filtered through a low-pass filter, applied in multiple step inputs in a 
sequential manner, and applied in multiple pulse inputs in a sequential 
manner. 
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